Damping of antiferromagnetic spin waves by valence fluctuations 
in the double layer perovskite YBaFe 2 5 
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Inelastic neutron scattering experiments show that spin dynamics in the charge ordered insulating 
ground state of the double-layer perovskite YBaFe2 0s is well described in terms of e g superexchange 
interactions. Above the Verwey transition at Ty = 308 K, t^g double exchange-type conduction 
within antiferromagnetic Fe02-BaO-FeC>2 double layers proceeds by an electron hopping process 
that requires a spin flip of the five-fold coordinated Fe ions, costing an energy of 5(J)S 2 « 0.1 eV. 
The hopping process disrupts near-neighbor spin correlations, leading to massive damping of zone- 
boundary spin waves. 

PACS numbers: 75.30.Ds,75.30.Et,72.25.-b,71.30.+h 



The extraordinary electrical transport properties of 
strongly correlated transition-metal oxides are often in- 
timately related to magnetic interactions. For example, 
in the doped perovskite manganites, electron spins in the 
itinerant e g band arc aligned parallel to the localized t 2g 
spins due to strong Hund's rule coupling. Consequently, 
conduction of spin polarized e g electrons results in fer- 
romagnetic (F) nearest-neighbor interactions, known as 
double exchange (DE) However, antiferromagnetic 
(AF) superexchange (SE) via localized t 2g electrons, co- 
operative Jahn-Teller distortions as well as onsite and in- 
tersite Coulomb repulsion may compete with DE leading 
to different insulating magnetic/orbital/charge-ordered 
phases @. EBaFeaOs (R = Nd - Ho and Y) dou- 
ble layer perovskites make up another class of compounds 
where DE and SE interactions are in competition, result- 
ing in either a charge-ordered insulator or valence-mixed 
conducting phase [J, [f|. However, in the i?BaFe205 
phases, DE occurs via the minority spin t 2g electrons 
in competition with AF SE dominated by the e g bands. 
Both SE and DE interactions may be modified in the 
charge/orbital-ordered ground state due to structural 
distortions and the fractional d xz orbital occupancy at 
the metal site [4(. Therefore, the study of spin excitations 
provides a window into the underlying charge dynamics. 
In this letter, we report on an inelastic neutron scatter- 
ing study of the magnetic excitations in YBaFe 2 05. The 
charge-ordered insulating ground state below the Verwey 
temperature Ty can be well understood in terms of nor- 
mal SE interactions. However, above Ty a strong cou- 
pling with valence fluctuations leads to massive damping 
of the zone-boundary spin waves. 

YBaFe205 is a fractional valent material with the 
Fe atoms possessing a nominal valence of 2.5+. The 
perovskite-based crystal structure of YBaFe205 con- 
sists of Fe02-BaO-FeC>2 double layers in which five- 
coordinated Fe sites form apex-shared square pyramids. 
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FIG. 1: (color online). Schematic representations of 
charge/magnetic order [J] in YBaFe2 0s. The upper and lower 
panels are projections onto the be and ab planes, respectively. 
Fe ions are shown as balls and oxygen square pyramids as 
triangles or squares. Arrows indicate magnetic moments, (a) 
For T < 308 K (AF C o), Fe 2+ (green online) and Fe 3+ (red 
online) chains run along the b axis and alternate along the a 
and c directions, (b) For 308 K < T < 430 K(AFvm) the Fe 
ions are valence mixed (blue online). 



The double layers are separated by an oxygen-vacant Y 
layer. At high temperatures, YBaFe20s is valence mixed 
and paramagnetic. On cooling, the compound orders an- 
tiferromagnetically at Tn = 425 K but remains valence 
mixed. At Ty — 308 K, a Verwey-type transition orders 
the Fe valences in real space and also modifies the AF 
structure. The magnetic/charge-ordering patterns 0] for 
T < Ty (AF CO ) and Ty < T < T N (AF VM ) are shown 
in FIG.[TJ Both AFco and AFvm structures are charac- 
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terized by strong AF coupling between Fe atoms within 
double layers due to e g SE interactions through nearly 
180° Fe-O-Fe bonds. The main difference between the 
two AF structures is the sign of the direct exchange inter- 
action between double layers (across the Y layer), which 
changes from F in AFvm to AF in AFco- However, 
this direct exchange is much weaker than AF SE within 
double layers due to the large Fe-Fe interlayer distance 
(~ 3.6 A 0]). This is shown using perturbation theory 
to estimate the ratio of direct exchange ( Jdiroct) and SE 
( Jaf) for a linear Fe-Fe bond @ 



■^direct tdd/U 

Jaf 



A 3 f&fc 



"Fc-O 



J10 

"Fc-Fc 



1%, (1) 



t^/A3 ft 4 u r,l da 

where t p d and tdd are the Fe-0 and Fe-Fe overlap inte- 
grals, U is the onsite Coulomb repulsion, and A is the 
charge transfer energy. These parameters are estimated 
as U = 8 eV and A = 3 eV 0] • Overlap integrals are esti- 
mated using Harrison's method [§] with dFc-Fc = 3.6 A, 
^Fc-o — 2.0 A and the coefficients of the radial over- 
lap integrals are r]ddc = 16.2 and r\^ va = 2.95 [8J]. 

To study the spin waves, we performed time-of-flight 
(TOF) powder inelastic neutron scattering measure- 
ments on approximately 50 g of YBaFe2 0s using the 
Pharos spectrometer at the Lujan Center, Los Alamos 
National Laboratory. The master sample of YBaFe2 0s 
was synthesized in several batches from citrate precursors 
as reported in Ref . Q • The combined batches were equi- 
librated at 400° C for 33 days in a sealed silica ampoule 
together with 0.03 g Zr foil as a getter, placed 12 cm from 
the sample and locally heated to 700° C. This was fol- 
lowed by a cool down at a rate of 0.1° C/min. The ceri- 
metrically determined oxygen content was 5.005(1) per 
formula at the top of this batch and 5.003(1) at the bot- 
tom. The sample quality was verified by X-ray powder 
diffraction, differential scanning calorimetry, and mag- 
netization measurements, which were in good agreement 
with previous results [4( . The powder sample was loaded 
in a flat aluminum can and mounted on the cold head 
of a closed-cycle He cryostat. For temperatures greater 
than 325 K, a displex with a high temperature stage was 
used. The sample was oriented at 45° to the incident 
neutron beam and inelastic neutron spectra were mea- 
sured at various temperatures between 6 K and 450 K, 
with an incident energy Ei = 120 meV. The TOF data 
were corrected for energy dependent detector efficiencies 
and instrumental background and reduced into energy 
transfer (E) and scattering angle (26) histograms. The 
resulting spectra contain both magnetic and phonon con- 
tributions. However, since the magnetic intensity disap- 
pears at high angles due to the magnetic form factor, it is 
possible to segregate the magnetic and phonon scattering 
by summing the data over different angle ranges. For the 
present experiment, we chose angle ranges 29 = 1 — 35° 
and 55 — 95° for low (magnetic) and high (phonon) an- 
gle data, respectively. The phonon contribution to the 
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FIG. 2: (color online), (a) Inelastic neutron scattering inten- 
sity, S(Q,E), for YBaFe 2 5 at T = 6 K. White lines are loci 
of constant angles in (Q, _E)-space that denote limits of the 
angle summation, 29 — 1° — 35°. (b) Calculation of S(Q,E) 
for YBaFe205 at T — 6 K using a Heisenberg model with 
J33 = 5.9 meV, J22 = 3.4 meV and J23 = 6.0 meV. (c) 
Comparison of angle-summed magnetic scattering data (open 
symbols) with Heisenberg model calculation (line). 



scattered intensity in the low angle data was removed by 
subtracting the appropriately scaled high angle data. 

The total scattered intensity for YBaFe205 at 6 K in 
the charge ordered AFco phase is shown as a function 
of Q (the magnitude of momentum transfer vector Q) 
and E in FIG. [2^a) . The main feature in the spectrum 
in the angle range 1° — 35° is a prominent band of scat- 
tering between 50 — 65 meV arising from zone-boundary 
spin waves. The band appears as two (or more) over- 
lapping peaks in the phonon subtracted magnetic inten- 
sity shown in FIG. [2c). Additional weak features below 
40 meV are due to imperfect subtraction of phonon bands 
near 20 meV and 40 meV, as well as dispersive features 
from both spin waves and phonons appearing as vertical 
streaks in FIG.^a). 

Measurements of polycrystalline-averaged spin excita- 
tion spectra of long-range magnetically ordered material 
result in neutron scattering intensities related to the spin- 
wave density-of-states (SWDOS). At low temperatures, 
the SWDOS of insulating and charge-ordered YBaFe20s 
may be understood by considering only SE interac- 
tions Q within double layers and is comparable to the 
G-type AF perovskite LaFe0 3 (see FIG.©. LaFe0 3 may 
be considered isotropic with a single nearest-neighbor ex- 
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FIG. 3: (color online). Magnetic scattering intensity 
for YBaFe 2 5 at 6 K and LaFe0 3 at 10 K. In- 
sets show idealized oxygen square pyramids(octahedra) in 
YBaFe2 05(LaFeOs) illustrating the broken magnetic bond in 
YBaFe2 05 in comparison to LaFe03. 



change constant Jo w 5 meV and spin S = 5/2 
Accordingly, the SWDOS of LaFeC>3 consists of a sin- 
gle sharp peak indicating a zone-boundary spin-wave en- 
ergy of 6JqS (~75 meV) due to the octahedral coordi- 
nation of the Fe sites. The average zone-boundary spin- 
wave energy for YBaFe205 (~60 meV) is very close to 
5J S. Therefore, the SWDOS in YBaFe 2 5 is compa- 
rable to LaFeC>3, after accounting for the one broken 
AF exchange bond (five nearest neighbors) present in the 
square-pyramidal coordination. 

While the comparison of the AFco spectrum to that 
of LaFeC>3 shows that the energy scale of the SWDOS 
in YBaFe205 is set by the average exchange (J) sa Jo, 
details such as the appearance of split peaks (see FIG. [2]) 
can only be explained by taking into account the varia- 
tions in the exchange between different Fe-Fe pairs. As 
illustrated in FIG.QJa), charge order in the AFco phase 
consists of Fe 2+ and Fe 3+ chains along the b axis, re- 
sulting in four unique nearest-neighbor SE paths, given 
by J33, J221 ^23 > an d J 23 5 where the subscripts indicate 
the valences of the Fe-Fe pairs and the superscripts de- 
note the direction of the magnetic bonds. It is reason- 
able to assume that the SE integral can be transferred 
from LaFe0 3 to the Fe 3+ -Fe 3+ pair in YBaFe 2 5 us- 
ing the approximate relationship j| 3 oc -^2 i 4 cos 2 6* [Io| . 
where t is the Fe 3+ -0 transfer integral (for 180° bond 
angles) which depends sensitively on the Fe-0 distance 
d as t oc d~ 7 / 2 [8] and 6 is the Fe-O-Fe bond angle. We 
have used the resulting estimate for J33 (5.9 meV) as a 
starting point for studies of model calculations using the 
Heisenberg Hamiltonian 
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where J^ is the exchange energy between spins Si and 
Sj and indicates that the sum is only over nearest 
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FIG. 4: (color online). Temperature dependence of the mag- 
netic scattering intensity for YBaFe2 0s between 6 K and 
450 K. Intensities at each temperature have been offset by 
a uniform amount for clarity. Spectra at 6 K (AFco) and 
316 K (AFvm) are compared in the inset. The lines are re- 
sults of a fit to a damped harmonic oscillator model(316 K 
data) and a Heisenberg model calculation (6 K data). 



neighbors. The Heisenberg Hamiltonian given in Eq. @ 
was used within linear spin-wave theory to obtain spin- 
wave energies and eigenvectors from which neutron in- 
tensities due to coherent scattering S(Q, E) were cal- 
culated. For a more detailed discussion of model spin- 
wave calculations, see Ref. ll|. Polycrystalline averaging 
of S(Q, E) was performed by Monte-Carlo integration 
over 10000 Q-vectors on each constant-Q sphere. The 
model calculations show that the two peak energies in 
YBaFe205 (at 58 and 62 meV) depend on 



E 2 
E 3 



2J22S2 



(2^23 

(2 J23 



J23) S3 

J23) ^2, 



(3) 



where Sj are 5/2 and 2 for Fe 3+ or Fe 2+ , respectively. 
These two energies are understood as the local excita- 
tion energy of the Fe 3+ (Fe 2+ ) spin inside the respective 
square pyramid. The combination J23 = (2Jf 3 + J 23 ) /3 
shifts the center-of-mass of the entire spin- wave band and 
J| 2 controls the splitting of the two peaks. Best agree- 
ment with both the splitting and the intensity ratio of the 
two peaks results in j| 3 w 5.9 meV, Jf 2 ~ 3.4 meV, and 
J23 ~ 6.0 meV for the charge-ordered state AFco- The 
calculated powder-averaged magnetic neutron intensities 
S(Q,E) and S(E), shown in FIG. [^b) and (c), respec- 
tively are in excellent agreement with the measured data. 

FIG. [4] shows the magnetic intensity as a function of 
temperature from 6 K to 450 K. Below T v = 308 K, the 
main spin- wave band is at 60 meV, as discussed above, al- 
though there is some broadening and softening in energy 
at 300 K, just below Ty. At 316 K, just above Ty in the 
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AFvm phase, a dramatic shift in the energy of the spin- 
wave peak down to 40 meV is accompanied by a huge 
broadening of the spin-wave band. Above Tn = 450 K, 
paramagnetic scattering is observed. The spectra in the 
AFco and AFvm phases are compared in the inset of 
FIG. [U The large decrease of zone-boundary spin- wave 
energy and the massive broadening in the valence mixed 
AFvm phase do not arise from normal SE interactions 
and can not be reproduced by a simple Heisenberg model. 
The broadening may arise from a distribution of exchange 
along different paths and/or from the presence of over- 
damped modes. In the AFvm phase, the valence-mixed 
Fe ions are electronically equivalent and occupy a single 
crystallographic site with Fe-O-Fe bond distances nearly 
identical in all directions. Thus, a decrease in the spin- 
wave bandwidth is expected in the AFvm phase in com- 
parison to AFco • 

The flipping of the sign of the direct exchange from 
AF (AFco) to F (AFvm) at Ty might appear to have a 
significant effect on the SWDOS. However, introduction 
of the direct exchange Jdirect to model spin-wave calcu- 
lations leads to a splitting of the SWDOS proportional 
to Jdirect rather than a decrease in energy of the entire 
spectrum. This confirms that Jdirect does not have a sig- 
nificant effect on the SWDOS, as was already expected 
from Eq. {Tj). Therefore, the broadening and decrease in 
energy of the zone-boundary spin waves suggest damp- 
ing arising from a strong coupling of spin excitations with 
valence-mixed charge carriers. 

The large changes in the spin-wave spectrum at Ty 
are explained by the appearance of new ferromagnetic 
interactions within the double layers. Valence- mixed 
tig bands give rise to ferromagnetic DE that competes 
with e g AF SE. In doped perovskite manganites, strong 
Hund's rule coupling of the itinerant e g electron to the 
localized ti g spin overwhelms the AF SE, leading to fer- 
romagnetism via DE In other words, metallic ferro- 
magnetism in the manganites is established in the limit 
Jde > Jse, where Jde and Jse represent DE and SE, 
respectively within a Heisenberg model. To estimate the 
strength of the DE within double layers, the spin-wave 
spectrum at 316 K in AFvm was fit using a damped har- 
monic oscillator model 

S ( E ) = \ _ e -E/k B T (£2 _ £2)2 + (Elf ' (4) 

where {l-e~ E / kB ' r )~ 1 is the Bose factor, F = 33(1) meV 
is the width and Eo = 45.8(4) meV is the position. The 
fit is shown as a red line in the inset of FIG. 01 We es- 
timate the average exchange energy from Eq = 5(J)S, 
resulting in (J) 4.07 meV; a 23% decrease from AFco- 
Given that the average exchange energy in the valence- 
mixed AFvm phase of YBaFe205 is (J) = Jse + Jde ~ 
4.07 meV and Jse ~ 5.3 meV in the charge-ordered 
AFco phase, we estimate that Jde ~ —1-2 meV and 
|Jde/Jse| ~ 0.23. In contrast to the manganites, con- 



duction in YBaFe205 proceeds via the minority spin t2 g 
electron, whose spin must be opposite that of the local 
spin on the half-filled Fe 3+ ion by the Pauli exclusion 
principle. Therefore, spin-polarized conduction within a 
double layer must overcome the AF structure established 
by the stronger, e g dominated SE. This can occur by 
flipping the spin on one Fe ion, momentarily creating an 
arrangement of five neighboring Fe spins that are ferro- 
magnetically aligned the central spin. It is then possible 
for real hopping of the ti g electron to occur within a dou- 
ble layer through a DE-type process. Flipping the spin 
at one Fe site costs an energy of 5(J)S 2 rs 0.1 eV. This 
scenario is consistent with conductivity measurements in 
GdBaFe205 [13| indicating that transport is an activated 
process in the valence mixed phase, with an activation en- 
ergy of about 0.1 eV. Such a hopping process will strongly 
disrupt near-neighbor spin correlations, leading to mas- 
sive damping of zone-boundary spin waves. From the 
perspective of the zone-boundary spin waves, DE implies 
that deviations of the relative angle between neighbor- 
ing spins makes electron hopping easier causing the 
decay of the zone-boundary spin wave into a charge and 
multi-magnon excitations. Given such a picture, a signif- 
icant change in the activation energy should be observed 
in resistivity in applied magnetic fields or at Tn- 
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